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Nitrogen is considered the most limiting nutrient for maize (Zea mays L), but there is limited understanding of the 
regulation of nitrogen-related genes during maize development. An Affymetrix 82K maize array was used to analyze the 
expression of <46 unique nitrogen uptake and assimilation probes in 50 maize tissues from seedling emergence to 31 d af- 
ter pollination. Four nitrogen-related expression clusters were identified in roots and shoots corresponding to, or overlap- 
ping, juvenile, adult, and reproductive phases of development. Quantitative real time PCR data was consistent with the ex- 
istence of these distinct expression clusters. Promoters corresponding to each cluster were screened for over-represented 
c/s-acting elements. The 8-bp distal motif of the /Arab/dops/s43-bp nitrogen response element (NRE) was over-represented 
in nitrogen-related maize gene promoters. This conserved motif, referred to here as NRE43-d8, was previously shown to be 
critical for nitrate-activated transcription of nitrate reductase (NIA1) and nitrite reductase (NIR1) by the NIN-LIKE PROTEIN 6 
(NLP6) in Arabidopsis. Here, NRE43-d8 was over-represented in the promoters of maize nitrate and ammonium transporter 
genes, specifically those that showed peak expression during early-stage vegetative development. This result predicts an 
expansion of the NRE-NLP6 regulon and suggests that it may have a developmental component in maize. We also report 
leaf expression of putative orthologs of nitrite transporters (NiTRI), a transporter not previously reported in maize. We 
conclude by discussing how each of the four transcriptional modules may be responsible for the different nitrogen uptake 
and assimilation requirements of leaves and roots at different stages of maize development. 



Introduction 

Plants have evolved multiple strategies to cope with wide vari- 
ation in the types and concentrations of soil nitrogen including 
nitrate and ammonium. 1 " 3 For example, at low external nitrate, 
plants employ high affinity transporters, while low affinity 
transporters are used when external nitrate is high. 1 Following 
uptake, soil-derived nitrogen is either transported directly to 
the shoot or first assimilated into amino acids for long-dis- 
tance transport and/or to provide nitrogen donor substrates for 
nitrogenous compounds. 4 To help adjust for changing nitrogen 
demand throughout plant development, 5 nitrogen incorporated 



in organic compounds can be scavenged from senescing tissues 
and re-converted into amino acids for transport through the 
phloem to growing tissues such as seeds. This process requires 
tissue specific expression of nitrogen assimilation enzymes and 
transporters. 4 To allow for such fine-tuned control, plants have 
evolved nitrogen-related gene families that perform similar func- 
tions but in different tissues and stages of development. 4,6 Such 
developmental control of nitrogen-related genes may involve 
distinct sets of transcriptional enhancers, but these are poorly 
characterized. In maize (Zea mays L.), very limited information 
has been reported on the transcriptional regulation of nitrogen 
uptake and assimilation genes across different tissues and stages 
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of development, which is therefore the objective of the current 
study. 

Nitrogen transporter genes have been most thoroughly char- 
acterized in Arabidopsis thaliana (L.). The Arabidopsis genome 
encodes at least 67 nitrate transporters, including 53 NRT1 
genes, 7 NRT2 genes, and 7 AtClc (chloride channel) genes. 6,7 
In maize, at least 4 genes encoding low affinity nitrate trans- 
porters (ZmNrtl.lA, ZmNrtl.lB, ZmNrtl.2, ZmNRT1.5A) have 
been reported, along with 4 genes encoding high affinity nitrate 
transporters (ZmNrt2.1, ZmNrt2.2, ZmNrt2.3, ZmNrt2.5). 2A ' n 
NRT2 requires interaction with co-transporter NAR2 (NRT3) 
proteins to be functionally active. 12 " 14 The maize genome encodes 
at least 2 NAR2-encoding genes [ZmNar2.1 (ZmNrt3.1), 
ZmNar2,2 (ZmNrt3.2)]. 2,4 In higher plants, a chloroplast-local- 
ized nitrite transporter (CsNitrl-L) was reported in cucumber, 
with a functional ortholog in Arabidopsis; 1 '' an NiTR ortholog 
has not been reported in maize. The genomes of higher plants 
have been reported to encode up to 14 ammonium transporter 
paralogs (AMT1 family), with the highest affinity transporters 
(AMTl;l; AMTl;3; AMT1;5) typically localized to root hairs 
and outer root cells; the lower affinity transporter(s) (AMTl;4) 
are located in the root endodermis. 4 In Arabidopsis, at least 5 
gene families have been reported to transport amino acids and 
peptides, some of which also transport inorganic nitrogen. 4 

Nitrogen assimilation is similarly complex. Higher plants 
directly assimilate ammonium into glutamine in plastids. 4 By 
contrast, nitrate is first converted to nitrite in the cytoplasm by 
nitrate reductase (NR), followed by reduction to ammonia in 
plastids by nitrite reductase (NiR). 4 Ammonium is fixed onto 
glutamate to form glutamine by glutamine synthetase (GS; Gin 
gene family), of which a plastidic isoform (GS2) and a cytosolic 
isoform (GS1) exist. A single gene in maize encodes GS2 {Gln2) 
whereas at least 5 genes encode GS1 (Glnl—1 to Glnl—5), which 
are differentially expressed during development. 16,17 Glutamine 
can subsequently react with 2-oxoglutarate to form 2 molecules 
of glutamate via glutamine 2-oxoglutarate amino transferase 
(GOGAT), also called glutamate synthase. 18 Glutamate can serve 
as an amino donor for other amino acids and nitrogen-requiring 
compounds, or act as an amine acceptor in the GS-GOGAT 
cycle to regenerate glutamine. Plants have 2 types of GOGAT 
enzymes, NADH-GOGAT and Fd-GOGAT, which use NADH 
and ferredoxin as electron donors, respectively. 18 Different 
GOGAT paralogs show constitutive or tissue-specific activity in 
plants, including in maize. 17,18 Ferredoxin-GOGAT is localized 
to leaf chloroplasts, while NADH-GOGAT is expressed in non- 
photosynthetic tissues including root plastids. 4 Following nitro- 
gen assimilation, glutamine, glutamate, and other amino acids, 
including asparagine and aspartate, are transported by vascular 
tissues to growing organs. 4 Nitrate and ammonium can also be 
stored in vacuoles before or after their long-distance transport. 4 

Transcriptional regulation of plant nitrogen-related genes is 
poorly understood, with very limited information available in 
monocots in particular. In dicots, including Arabidopsis, a 43 
bp pseudo-palindromic nitrogen response element (NRE) was 
recently identified in the promoter of a nitrite reductase (NiR) 
gene and shown to be necessary and sufficient for nitrate-activated 



gene transcription. 19 Subsequent research showed that a similar 
43 bp NRE is present in the promoters of NiR genes from rice, 
maize, and sorghum, and that 4 copies of the maize NRE are 
necessary and sufficient for nitrate induction of a GUS reporter 
construct in Arabidopsis. 20 The authors reported that they did 
not observe the NRE in genes other than nitrate reductase 
(NIA1) and nitrite reductase {NIR1) in Arabidopsis, 21 consistent 
with a parallel study. 22 NRE-like sequences were also identified 
downstream of the transcriptional terminator of the Arabidopsis 
nitrate reductase NIA1 gene. 23 A different research group, how- 
ever, reported the existence of 3 non-NRE motifs that were 
necessary and sufficient for nitrate induction of the Arabidopsis 
NIA1 promoter, specifically: a myb motif (C-G/C-GTT-G/A), 
ALFIN1 (core C/A-CAC), and an E-box (CAAGTG). 24 

Plant development above ground has been divided into dis- 
tinct developmental phases based on characteristics including 
leaf shape, surface wax, the presence of trichomes, and under- 
lying genetic networks. 25,26 In maize, the juvenile vegetative 
phase (V) is from seedling emergence (Ve) to the formation of 
leaves 4 to 6 (V2 stage), which is followed by the adult vegetative 
phase (here > V5). Flowering initiates the reproductive phase 
(R), which includes seed formation, the latter divided into stages 
based on the number of days after pollination (DAP). Much less 
information has been reported on whether these phases have any 
importance in categorizing root development. 

The objective of this study was to characterize the global 
expression pattern of nitrogen transport and assimilation genes 
during the different phases of maize development, and to iden- 
tify associated aV-element motifs over-represented in the pro- 
moters corresponding to each expression cluster. We analyzed 
data from a developmental expression series which we recently 
reported for maize. 27 Specifically, the expression of 46 putatively 
unique nitrogen uptake and assimilation probes were analyzed 
in 50 maize tissues from seedling emergence to 31 d after pol- 
lination using an 82K customized maize B73 inbred microarray, 
with a focus on root and leaf expression. The expression of a 
subset of these probes was subsequently validated using quan- 
titative real time PCR (qRT-PCR). Promoters associated with 
distinct expression clusters were subsequently used to screen for 
over-represented cw-acting elements using a new de novo motif 
discovery program customized and validated for maize called 
Promzea (Promoters of Zed). 28 

Results 

Tissue and developmental stage selective expression 

As we recently described, 27 RNA was collected from 50 tis- 
sues from seedling emergence to 31 d after pollination (DAP; 
Table SI) and hybridized onto a customized maize Affymetrix 
82K Unigene array. From all the probes present on the array, 
and after substantial trimming (Table S2; see below for details), 
a total of 65 nitrogen uptake and assimilation array probes 
were identified (Table 1) and analyzed, of which 46 probes 
were predicted to be unique. The expression patterns for these 
nitrogen-related probes from all 50-tissue samples are sum- 
marized (Fig. SI). Expression of 64 of the 65 probes showed 
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significant differential expression in at least 1 of the tissue-stages 
of development, with adjusted p-values ranging from 4.6le" 04 
to 4.51e" 36 , well below the cut-off of 0.05 (data not shown). 
Only 1 probe, ZmN-25, was not differentially expressed; this 
probe corresponded to an oligopeptide and nitrate transporter 
(GRMZM2G127134) in expression Cluster 2. 

In reproductive tissues, consistently high nitrogen-related 
expression was observed in anthers (Fig. SI). Interestingly, 
probes corresponding to nitrite transporters (NiTRl), 15 a gene 
class not previously reported in maize, were expressed at differ- 
ent stages of maize development (Fig. 1; Fig. SI). In vegetative 
tissues, 4 nitrogen-related gene expression clusters were detected 
(Fig. 1; Table S3). Cluster 1 probes showed peak expression in 
juvenile stage roots (Ve-V2), with several probes also peaking 
in vegetative stage leaves (Ve-V5), suggestive of a whole plant 
juvenility cluster. Cluster 1 probes matched genes encoding a 
high affinity transporter NRT2.1, the low affinity transporter 
NRTl.l, and ammonium transporters (Table 1). Cluster 2 
probes were leaf-selective and more highly expressed in vegeta- 
tive stage leaves (juvenile and adult: Ve-V5) than post-flowering 
leaves (R1-R31) (Fig. 1; Table S3). Cluster 2 probes matched 
low affinity nitrate transporters including NRT1.5, nitrite 
transporter NiTRl, and nitrate reductases (NR1, NR2) (Table 
1). Cluster 3 probes were root selective and were expressed 
throughout development (juvenile to post-flowering) (Fig. 1; 
Table S3). Cluster 3 probes matched genes encoding glutamine 
synthetases (cytosolic isoform GS1) including Gln4IGlnl—4 
and/or Gln5IGlnl—5 as well as several glutamate synthases 
(GOGAT, Table 1). Also included in Cluster 3 were ammo- 
nium transporters, the low affinity transporter NRTl.l, and 
the companion protein (NAR2.1) of the high affinity nitrate 
transporter complex. Finally, Cluster 4 probes were leaf selec- 
tive and showed the most consistent expression in older leaves at 
later stages of development (adult to post flowering: V5-R31), 
with 3 probes showing peak expression in reproductive stage 
roots (Fig. 1; Table S3). Several Cluster 4 probes matched genes 
encoding nitrate reductase (NR1) and nitrate transporters, as 
well as 1 NAR2 paralog and the nitrite transporter NiTRl 
(Table 1). 

Validation of microarray expression 

Using corresponding gene sequences from the maize genome, 
expression of 8 nitrogen-related probes was verified by quan- 
titative real time PCR (qRT-PCR) for 6 of the developmental 
tissue-stages (Fig. 2; Fig. S2; Tables S4 and S5). The results 
showed that 75% of the qRT-PCR data sets (6 out of 8) cor- 
related with the microarray results, having Pearson coefficients 
of correlation (r) ranging from 0.8 to 0.97 (Fig. S2; Table S5). 
The Pearson coefficients of correlation were significant for 5 of 
these 6 genes (/"value < 0.05) while the sixth showed borderline 
significance (/"value = 0.056); 2 were not correlated {Nrtl.l and 
AW) (Fig. S2; Table S5). Specifically, the high affinity nitrate 
transporter Nrt2.1 showed peak expression in VI seminal roots, 
and a minor peak in V5 leaves (Fig. 2A), consistent with belong- 
ing to juvenility Cluster 1 (Fig. 1). Though not statistically 
correlated with Cluster 1, expression of low affinity transporter 
Nrtl.l also showed peak expression in VI seminal roots (Fig. 



2B). Expression of nitrate reductase Nr2 showed leaf-selective 
expression (note expanded y-axis values), peaking in VI leaves 
(Fig. 2C), consistent with belonging to Cluster 2. Though pre- 
dicted by microarray data to belong to Cluster 2, qRT-PCR 
expression of Nrl appeared to be more consistent with juvenil- 
ity Cluster 1, with expression peaks in both juvenile roots and 
leaves, and little expression in reproductive-staged leaves and 
roots (Fig. 2D). Nar2.1 showed peak expression in VI seminal 
roots, and higher expression in the 3 root tissues compared with 
the 3 leaf tissues examined (Fig. 2E), consistent with belonging 
to Cluster 3. Cluster 3 was also exemplified by NADH glutamate 
synthase (NADH-GOGAT), which showed very high expression 
in VI and V5 seminal roots, but low expression at all leaf stages 
(Fig. 2F). Nar2.2 showed peak expression in reproductive stage 
(Rl) leaves but low expression in Rl roots (Fig. 2G), consistent 
with belonging to reproductive stage Cluster 4. A high affin- 
ity nitrate transporter (GRMZM2G455124) also showed peak 
expression in Rl stage leaves, but also a second peak in Rl nodal 
roots (Fig. 2H), consistent with a reproductive stage expression 
cluster. 

More generally, the qRT-PCR expression data was consistent 
with maize having distinct nitrogen gene expression clusters 
during vegetative development; this became apparent when pairs 
of genes were compared (Fig. 2). For example, Nr2 and NADH 
glutamate synthase showed opposite expression, peaking in juve- 
nile leaves vs. juvenile roots, respectively (Fig. 2C, D). Similarly, 
expression of a high affinity nitrate transporter (Maize GDB ID: 
GRMZM2G455124), which peaked in reproductive-stage veg- 
etative tissues (Fig. 2H), was opposite to Nrl which peaked in 
juvenile stage tissues (Fig. 2D). Finally, Nar2.1, which showed 
peak expression in VI seminal roots and exhibited low expres- 
sion in reproductive stage (Rl) leaves (Fig. 2E), was strikingly 
opposite to Nar2.2 which showed peak expression in Rl leaves 
(Fig. 2G). 

Over-represented cw-acting promoter motifs 

Four over-represented motifs were present in the promoters of 
the nitrogen-related genes expressed in Cluster 1 (early root/leaf 
cluster), of which 3 were similar to previously identified oV-act- 
ing elements (Table 2). Interestingly, 1 Cluster 1 motif was very 
similar to the conserved 8 bp distal element of the 43 bp pseudo- 
palindromic nitrogen response element [(t/c)GACcCTT] 
identified in dicot and monocot nitrite reductase (NiR) genes, 
including maize, and shown to be necessary and sufficient for 
nitrate-activated gene transcription. 19,20 Within Cluster 1, this 
8 bp motif (consensus CGACCCTT), which we refer to as 
NRE43-d8 (NRE 43 bp, distal 8 bp, see Discussion), was over- 
represented in the promoters of genes encoding nitrate trans- 
porters (NRT1, NRT2.1) and ammonium transporters (Table 
3). The second Cluster 1 element was similar to the G-box 10 
and E-Box elements (Table 2). The core of the G-box 10 ele- 
ment is identical to the E-box element (CANNTG) which was 
recently shown to be critical for nitrate induction of the NIA1 
(nitrate reductase) promoter in Arabidopsis. 24 The third Cluster 
1 element, NRG1 (Negative Regulator of Glucose-repressed 
Genes 1) (Table 2), was previously identified in yeast upstream 
of genes involved in regulating cell sugar status. 29,30 For example, 
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ZmN-6 ammonium transporter 

ZmN-1 1 ammonium transporter (Osamd .1 ) 

ZmN-1 high affinity nitrate transporter (nrt2.1) 

ZmN-1 3 NADH-dependent glutamate synthase 

ZmN-4 nitrate and chloride transporter 

ZmN-5 nitrate and chloride transporter 

ZmN-9 nitrate transporter (nrt1 . 1 ) 

ZmN-2 nitrate transporter (nrt1 . 1 ) 

ZmN-3 nitrate transporter (nrt1.1) 

ZmN-12 nitrate transporter (nrt1.1) 

ZmN-7 symbiotic ammonium transporter 

ZmN-1 0 symbiotic ammonium transporter 

ZmN-8 symbiotic ammonium transporter 

ZmN-28 ferredoxin-dependent glutamate synthase 
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ZmN-1 9 nitrate and chloride transporter 
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ZmN-22 nitrate transporter (nrll -5) 

ZmN-23 nitrite transporter (NiTR1 gene) 

ZmN-25 oligopeptide and nitrate transporter 

ZmN-29 amino acid transport protein 

ZmN-39 ammonium transporter 

ZmN-32 ammonium transporter 

ZmN-37 ammonium transporter 

ZmN-34 ammonium transporter (OsAMT! . 1 ; Osamt5.2) 

ZmN-40 gln4; gln5 

ZmN-41 gln4; gln5 

ZmN-30 high affinity nitrate transporter(nar2.1 ) 

ZmN-42 NADH-dependent glutamate synthase 

ZmN-44 NADH-dependent glutamate synthase 

ZmN-43 NADH-dependent glutamate synthase 

ZmN-45 NADH-dependent glutamate synthase 

ZmN-46 NADH-dependent glutamate synthase 

ZmN-36 nitrate reductase (nr2) 

ZmN-35 nitrate transporter 

ZmN-31 nitrate transporter (nrtl.1) 

ZmN-33 nitrate transporter (nrt1 .1 ) 

ZmN-38 oligopeptide and nitrate transporter 

ZmN-58 ammonium transporter 

ZmN-61 ammonium transporter (OsAMT3.3) 

ZmN-55 high affinity nitrate transporter 

ZmN-56 high affinity nitrate transporter (nar2.1 ; nar2.2) 

ZmN-65 NADH-dependent glutamate synthase 

ZmN-48 nitrate reductase (NADH; nrl) 

ZmN-49 nitrate reductase (NADH; nrl) 

ZmN-52 nitrate reductase (NADH; nrl) 

ZmN-60 nitrate reductase (NADH; nrl) 

ZmN-64 nitrate reductase (NADH; nM) 

ZmN-47 nitrate transporter 

ZmN-54 nitrate transporter 

ZmN-57 nitrate transporter 

ZmN-63 nitrate transporter 

ZmN-51 nitrate transporter 

ZmN-50 nitrate transporter (AtntM) 

ZmN-53 nitrate transporter (Glycine max nrt1 .2) 

ZmN-59 nitrite transporter (NiTR1 gene) 

ZmN-62 nitrite transporter (NiTR1 gene) 
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Figure 1. Heatmap showing relative expression of nitrogen-related probes in root and leaf tissues at different developmental stages. Each row rep- 
resents an array probe set. The relative expression level of each probe set ranges from low expression (yellow) to high expression (dark blue). V1-V5 
stages refer to the number of visible stem-leaf nodes. Tissue sampling was as follows: Ve leaf: coleoptile tissue; VI leaf: leaves 1 and 2; V2 leaf: actively 
growing leaf 4; V5 leaf: actively growing leaf 8, 15 cm from the tip; R1-R31 leaf: second leaf above top ear, 15 cm from the tip; VI -V5 sroot: seminal root 
from vegetative stages; V2-V5 nroot: nodal root (crown root) from vegetative stages; R1-R24 nroot: nodal root (crown root) from reproductive stages. Ab- 
breviations: V, vegetative pre-flowering stage; R, reproductive post-flowering stage; nDAP, days after pollination. For additional details, refer to Table 1; 
Fig. SI; Tables S1 and S3. 



the NRGl-like motif was over-represented in the promoters of With respect to the additional expression clusters, there were 
Nrtl.l and Nrt2.1 (data not shown). two over-represented motifs in the promoters of Cluster 2 (early 
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Table 1. Microarray probes used in this study 



Expression 
cluster 


Probe 
number 


GenBank ID 


Maize GDB ID 


Corresponding gene annotation 


4 


ZmN-65 


BAA35120.1 


GRMZM2G077054; GRMZM2G085078; GRMZM2G375064 


NADH-dependent glutamate synthase 


4 


ZmN-64 


AAD38068.1 


GRMZM2G076723; GRMZM2G443637 


Nitrate reductase (NADH; nrl) 


4 


ZmN-63* 


BAB89087.1 


GRMZM2G064091 


Nitrate transporter 


4 


ZmN-62 


CAA93316.2 


GRMZM5G827496 


Nitrite transporter (NiTR1 gene) 


4 


ZmN-61 


AAN06953.1 


GRMZM2G043193; GRMZM2G080045; GRMZM2G338809 


Ammonium transporter (OsAMT3.3) 


4 


ZmN-60 


AAD38068.1 


GRMZM2G076723; GRMZM2G443637 


Nitrate reductase (NADH; nrl) 


4 


ZmN-59 


CAA93316.2 


GRMZM5G827496 


Nitrite transporter (NiTR1 gene) 


4 


ZmN-58 


BAC65231.1 


GRMZM2G080045 


Ammonium transporter 


4 


ZmN-57* 


BAB89087.1 


GRMZM2G064091 


Nitrate transporter 


4 


ZmN-56 


AAM12786.1 


GRMZM2G179294; GRMZM2G163494 


High affinity nitrate transporter (nar2.1; nar2.2) 


4 


ZmN-55 


AAF78499.1 


GRMZM2G455124 


High affinity nitrate transporter 


4 


ZmN-54* 


BAB89087.1 


GRMZM2G064091 


Nitrate transporter 


4 


ZmN-53 


BAB1 9757.1 


GRMZM2G122251 


Nitrate transporter (Glycine max nrt1.2) 


4 


ZmN-52A 


AAD38068.1 


GRMZM2G076723; GRMZM2G443637 


Nitrate reductase (NADH; nrl) 


4 


ZmN-51 


CAD41 034.1 


GRMZM2G176253 


Nitrate transporter 


4 


ZmN-50 


AAG52554.1 


GRMZM2G137421 


Nitrate transporter (Atntll) 


4 


ZmN-49^ 


AAD38068.1 


GRMZM2G076723; GRMZM2G443637 


Nitrate reductase (NADH; nrl) 


4 


ZmN-48A 


AAD38068.1 


GRMZM2G076723; GRMZM2G443637 


Nitrate reductase (NADH; nrl) 


4 


ZmN-47 


AAO00921.1 


GRMZM2G137421 


Nitrate transporter 












3 


ZmN-46* 


BAC01257.1 


GRMZM2G077054; GRMZM2G085078; GRMZM2G375064 


NADH-dependent glutamate synthase 


3 


ZmN-45* 


BAC01 257.1 


GRMZM2G077054; GRMZM2G085078; GRMZM2G375064 


NADH-dependent glutamate synthase 


3 


ZmN-44 


BAA35120.1 


GRMZM2G077054; GRMZM2G085078; GRMZM2G375064 


NADH-dependent glutamate synthase 


3 


ZmN-43^ 


BAC01257.1 


GRMZM2G077054; GRMZM2G085078; GRMZM2G375064 


NADH-dependent glutamate synthase 


3 


ZmN-42A 


BAA35120.1 


GRMZM2G077054; GRMZM2G085078; GRMZM2G375064 


NADH-dependent glutamate synthase 


3 


ZmN-41# 


P38561 


GRMZM5G872068; GRMZM2G036464 


Glutamine synthase gln4; gln5 


3 


ZmN-40# 


P38561 


GRMZM5G872068; GRMZM2G036464 


Glutamine synthase gln4; gln5 


3 


ZmN-39 


AAL05614.1 


GRMZM2G028736 


Ammonium transporter 


3 


ZmN-38¥ 


AAD39317.1 


GRMZM2G127134 


Oligopeptide and nitrate transporter 


3 


ZmN- 

37** 


BAC65231.1 


GRMZM2G080045 


Ammonium transporter 


3 


ZmN-36 


P27968 


GRMZM5G878558 


Nitrate reductase (nr2) 


3 


ZmN-35¥ 


AAD39317.1 


GRMZM2G127134 


Nitrate transporter 


3 


ZmN-34 


AAL05612.1 


GRMZM2G1 18950; GRMZM2G473697; GRMZM2G701289 


Ammonium transporter (OsAMT1.1; Osamt5.2) 


3 


ZmN-33t 


BAC56914.1 


GRMZM2G086496 


Nitrate transporter (nrt1 .1) 


3 


ZmN- 

32** 


BAC65231.1 


GRMZM2G080045 


Ammonium transporter 


3 


ZmN-31t 


BAC56914.1 


GRMZM2G086496 


Nitrate transporter (nrtl .1) 


3 


ZmN-30 


AAM12786.1 


GRMZM2G179294 


High affinity nitrate transporter (nar2.1) 


3 


ZmN-29 


AAB82307.1 


GRMZM2G1 73967 


Amino acid transport protein 












2 


ZmN-28 


MZEFEGLU 


GRMZM2G036609 


Ferredoxin-dependent glutamate synthase 


2 


ZmN-27* 


P27968 


GRMZM5G878558 


Nitrate reductase (nr2) 
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Expression 
cluster 


Probe 
number 


GenBank ID 


Maize GDB ID 


Corresponding gene annotation 


2 


ZmN-26 


P27967 


GRMZM2G568636 


Nitrate reductase (NADH) 


2 


ZmN-25 


AAD39317.1 


GRMZM2G127134 


Oligopeptide and nitrate transporter 


2 


ZmN-24A 


AAG52554.1 


GRMZM2G137421 


Nitrate transporter (Atntll) 


2 


ZmN-23 


CAA93316.2 


GRMZM5G827496 


Nitrite transporter (NiTR1 gene) 


2 


ZmN-22 


CAD41 141.1 


GRMZM2G138731 


Nitrate transporter (nrtl-5) 


2 


ZmN-21* 


P27968 


GRMZM5G878558 


Nitrate reductase (nr2) 


2 


ZmN-20A 


AAG52554.1 


GRMZM2G137421 


Nitrate transporter (Atntll) 


2 


ZmN-19t 


NoDescription 


NA 


Nitrate and chloride transporter 


2 


ZmN-18 


P49102 


GRMZM2G568636 


Nitrate reductase (NADH) 


2 


ZmN-17 


AAD38068.1 


GRMZM2G076723; GRMZM2G443637 


Nitrate reductase (NADH; nrl) 


2 


ZmN-16 


BAB89087.1 


GRMZM2G064091 


Nitrate transporter 


2 


ZmN-15t 


AAL31925.1 


GRMZM2G453320 


Nitrate and chloride transporter 


2 


ZmN-14t 


AAL31925.1 


GRMZM2G453320 


Nitrate and chloride transporter 












1 


ZmN-13 


BAA35120.1 


GRMZM2G077054; GRMZM2G085078; GRMZM2G375064 


NADH-dependent glutamate synthase 


1 


ZmN-12 


Q05085 


GRMZM2G086496 


Nitrate transporter (nrt1 .1) 


1 


ZmN-11 


AAL05612.1 


AC208641 .3_FG005; AC208641 .3_FG008; 
GRMZM2G028736; GRMZM2G175140; GRMZM2G335218 


Ammonium transporter (Osamtl.1) 


1 


ZmN-10* 


No 

Description 


NA 


Symbiotic ammonium transporter 


1 


ZmN-9A 


BAC56913.1 


GRMZM2G086496 


Nitrate transporter (nrtl .1) 


1 


ZmN-8* 


AAC32828.1 


GRMZM2G082343 


Symbiotic ammonium transporter 


1 


ZmN-7* 


AAC32828.1 


AC211704.3_FG003; GRMZM2G082343 


Symbiotic ammonium transporter 


1 


ZmN-6 


BAC65231.1 


GRMZM2G080045 


Ammonium transporter 


1 


ZmN-5t 


AAL31925.1 


GRMZM2G453320 


Nitrate and chloride transporter 




ZmN-4t 


AAL31925.1 


GRMZM2G453320 


Nitrate and chloride transporter 




ZmN-3A 


BAC56914.1 


GRMZM2G086496 


Nitrate transporter (nrtl .1) 




ZmN-2A 


BAC56914.1 


GRMZM2G086496 


Nitrate transporter (nrtl .1) 




ZmN-1 


AY129953 


GRMZM2G010280 


High affinity nitrate transporter (nrt2.1) 



Note: Within each cluster, probe numbers that share a common symbol are predicted to correspond to the same gene based on the sequence homology 
of the corresponding genes and their similar expression patterns. 



leaf cluster): LTRE (low temperature responsive element) 31 and 
YAP1 (Table 2). 32 ' 33 The core of the LTRE element, also called 
C/DRE, mediates ABA-independent responses to cold, ampli- 
fied by light via phytochrome signaling. 31 The YAP1 element 
was first identified as the binding site of the yeast YAP1 protein, 
a conserved protein across eukaryotes including Arabidopsis? 4 
YAP1 was shown to regulate the yeast cell cycle in response 
to oxidative stress, but also able to regulate cell elongation. 32 ' 33 
Over-represented in Cluster 3 promoters (root cluster) was a 
motif similar to the binding site of the Forkhead (FHL1) tran- 
scription factor from yeast (Table 2), a repressor of ribosome- 
encoding genes during glucose starvation. 35 Here, the FHL1 
motif was over-represented in genes encoding nitrate and ammo- 
nium transporters, nitrate reductase, and glutamine synthases 
(data not shown). Finally, a previously unidentified motif was 



found in the promoters of Cluster 4 genes (late leaf cluster) with 
the consensus NANGAG (Table 2). 

Discussion 

Nitrogen-related genes have distinct patterns of expression 
during maize vegetative development 

Recent transcriptome studies have highlighted global gene 
expression patterns during maize development, examining 
gene expression along a maturation gradient of a juvenile leaf; 5 
in leaves during the transition from juvenile to adult vegetative 
stages; 36 and in different organs and tissues across developmental 
stages. 37 ' 27 The focus of this study was to examine the expres- 
sion of nitrogen transporter and assimilation genes across maize 
development, particularly in roots and leaves. Four vegetative 
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gene expression clusters were identified (Fig. 1; Table S3), 
specifically an early stage-selective cluster (Cluster 1, roots 
and leaves), an early leaf-selective cluster (Cluster 2, juvenile 
to adult leaves), a root-selective cluster (Cluster 3, juvenile to 
reproductive roots), and a late leaf-selective cluster (Cluster 4, 
adult to reproductive leaves), which also included some late- 
stage root expression. We observed that, within a gene family, 
distinct subsets of nitrogen-related paralogs were preferentially 
expressed in leaves compared with roots (Fig. 1; Fig. 2; Table 
S3). For example, high affinity co-transporter Nar2.1, showed 
peak expression in roots (Fig. 2E), while paralog Nar2.2 
showed peak expression in leaves (Fig. 2G). Furthermore, dis- 
tinct paralogs were expressed at different developmental phases 
corresponding to previously defined juvenile, adult and repro- 
ductive phases of shoot development in maize. 25,26 For example, 
qRT-PCR data showed that a high affinity nitrate transporter 
(Nrt2.1) was expressed in juvenile-stage roots and leaves (Fig. 
2A) compared with a different high affinity nitrate transporter 
(GRMZM2G455124) which was expressed in reproductive- 
stage roots and leaves (Fig. 2H). It is noteworthy that nitro- 
gen-related gene expression in roots also corresponded to these 
aboveground developmental shifts (Fig. 1; Table S3). The latter 
result may not be surprising, however, given that earlier stud- 
ies have demonstrated that nitrogen stress responses are highly 
coordinated between the shoot and root systems of maize. 2,38,39 

As a note added in proof, our root expression results are 
consistent with a recently published study by Garnett et al. 11 
which analyzed patterns of expression for nitrate uptake trans- 
porter paralogs in roots during maize development. Garnett 
et al. showed that low affinity {Nrtl family) and high affin- 
ity (Nrt2, Nar2/Nrt3 family) nitrate transporters change in 
expression during maize development: peak abundance of 
nitrate transporter transcripts generally occurs in juvenile and 
reproductive stage roots. Furthermore, Garnett et al. showed 
that within each transporter gene family, different paralogs 
exhibit different expression patterns: for example ZmNRTl.lA 
peaks early in development in roots, whereas ZmNRT1.2 peaks 
in roots during reproductive stages." Direct comparisons are 
difficult to make between our study and Garnett et al., how- 
ever, because the latter study employed inbred Gaspe Flint 
which has a much shorter life cycle and dramatically smaller 
shoot compared with the Syngeta inbred used here; the two 
studies also employed different nitrogen and growth condi- 
tions (hydroponics vs. coarse clay). 

The developmental complexity of nitrogen-related gene 
expression observed in our study is also consistent with data 
emerging from other species. For example, in Arabidopsis, 
Nrtl gene family paralogs are differentially expressed in roots, 
stems, leaves and reproductive tissues. 6,40 Also in Arabidopsis, 
nitrogen-related genes including nitrite reductase (NIR1), 
nitrate reductase (NIA1), and a high affinity nitrate trans- 
porter (Nrt2A), showed differential expression in young roots 
vs. young leaves during nitrogen starvation. 41 

High expression in anthers 

In terms of reproductive development, a particularly interest- 
ing finding from this study was the high level of gene expression 
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Figure 2. Quantitative real time PCR (qPCR) analysis of select genes in 
order to verify: microarray profile results, and the existence of distinct 
clusters of gene expression of nitrogen-related genes in maize. The labels 
on the left indicate the original microarray cluster predictions. The y-axis 
shows the relative level of gene expression; please note the different scales 
used. Shown are results for: (A) high affinity nitrate transporter Nrt2.1, (B) 
low affinity nitrate transporter Nrf?. 1, (C,D) nitrate reductases (Nr2, Nr1), (E) 
co-transporter for high affinity nitrate transporter {Nar2.1), (F) NADH-GO- 
GAT, (G) co-transporter for high affinity nitrate transporter (Nar2.2), and (G) 
a high affinity nitrate transporter (GRMZM2G455124). The error bar repre- 
sents the standard error of the mean of 3 biological replicates. Please refer 
to Table S4 for gene information and PCR primers used, and to Fig. S2 and 
Table S5 for comparisons between qPCR and microarray results. 



of many nitrogen-related probes in anthers (Fig. SI). In fact, 
> 14/65 probes showed peak expression in anthers (VT stage), 
including probes that showed root or shoot-selective expression 
during vegetative development. Given that the normalization 
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Table 2. Putative c/s-acting promoter motifs over-represented in the nitrogen-related genes expressed during maize development at different stages/ 
tissues 



Cluster 


Motif Logo 


Consensus 


Preferential position 


Similarity 


E-value 


Alignment 






CGACCNTT 


-50 +1 


NRF rnrp 


2.26E-09 


CGACCNTT 
YGACCCTT 




* lcCA TGc 


CCACGTGC 


-300 -250 


G-BOX 10 
reverse 


8.78E-11 


-GCANCNTGG- 
GGCACGTGGC 


1 


"\caC cC Ax 

o 1 _ „ 5. „ _ . Sr 2T i ^ 


CACGCCRCAG 


-250 -200 


/ 


/ 


/ 






NNANGCSGCW 


-50 +1 


NRG1 reverse 


6.23E-07 


AANGGTCG 
-AGGTCC 


2 




AGTCGG 


-500 -450/ -350 -300 


LTRE 


6.88E-07 


AGTCGG 
-GTCGG 


■•1 _ _ ft c 


TAGYCRGC 


-250 -200 


YAP1 


6.73E-07 


— TAGYCRGC- 
NTTAGTMAGCN 


3 




TRTCCGTACG 


-350 -300 


FHL1 


1.19E-06 


TRTCCGTACG- 
-AYOCGTACAT 


4 




NANGAG 


-500 -450 


/ 


/ 


/ 



Gene expression clusters correspond to Figure I.The most common position(s) of the motif in the -500 bp promoter region is indicated (Preferential 
position). Matches to previously identified motifs in promoter databases are identified (Similarity) along with the similarity score (E-value) and alignment 
to the match (Alignment). Only the most over-represented motifs are shown. Abbreviations: n = any nucleotide, Y = C/T, R = A/G, S = G/C, W = A/T. 



method used, named RMA, 42 standardizes the gene expression 
ranges between arrays (e.g., relative gene expression from anther 
RNA has the same range than leaf RNA after normalization), 
this result does not appear to be an artifact. It is also noteworthy 
that probes corresponding to nitrate transporters were expressed 
in late stage tassels (Fig. SI). 
Nitrite transporter expression 

As noted in the above section, transcripts of putative maize 
nitrite transporter (NiTR) orthologs were detected in the 2 
leaf-selective expression clusters (Clusters 2 and 4; Fig. 1) with 
additional strong expression in the husk leaves surrounding the 
cob (Fig. SI; Table S3) but were excluded from Clusters 1 and 3 
which included root expressed genes. This data is consistent with 
nitrite transporters being localized to chloroplasts. 6 Here we note 
NiTRs in a separate section because this gene class has not pre- 
viously been reported in maize. A gene encoding NiTR (Narl) 
was initially reported in Chlamydomonas chloroplasts. 43 In higher 
plants, a NiTR gene was first reported in cucumber (CsNitrl-L), 
along with a functional ortholog tested inArabidopsis (Atlg68570) 
where a knockout mutation showed a 5-fold increase in nitrite 
accumulation in leaves. 15 The cucumber NiTR protein was local- 
ized to the inner envelope membrane of chloroplasts, where it was 
hypothesized to load nitrite from the cytoplasm into the stroma 
of the chloroplast during nitrate assimilation." Functional data 
will be needed to authenticate the genes corresponding to the 
putative maize NiTR(s) probes reported in this study. 



Study limitations including discrepancies between microar- 
ray data and qRT-PCR data 

In this study, the reader is cautioned not to interpret each 
probe as representing a unique locus. The microarray was 
designed using EST sequences prior to the major sequence 
release of the maize genome. As a result, more than 1 nitrogen- 
related probe could match the same gene; this scenario was pos- 
sible for a total of 49 probe sets (75%; 49/65). These related 
probe sets could have originated from different portions of the 
same gene including alternative splice products (e.g., probes 
ZmN-2 and ZmN-3 which have the same array expression pat- 
tern, Table S3). Alternatively, the related probe sets may have 
originated from closely related paralogs due to the extensive 
number of duplications in the allotetraploid maize genome (e.g., 
nitrate transporter probes ZmN-31 and ZmN-33 have a differ- 
ent array expression pattern than ZmN-2 and ZmN-3) (Fig. 1; 
Table S3). The latter hypothesis could explain the low corre- 
lation observed for 2 microarray probes when compared with 
qRT-PCR data, specifically the probes corresponding to Nrtl.l 
in Cluster 2 (probe ZmN-2) and Nrl in Cluster 3 (Fig. S2). The 
microarray vs. qRT-PCR data might in fact represent duplicated 
paralogs (Table S4). Recent gene duplication phenomena (e.g., 
tandem duplications) were suggested by the observation that 
18 nitrogen-related probes (27%; 18/65) matched more than 1 
maize locus (Table 1). Based on sequence identities and expres- 
sion patterns, we estimate that the maximum number of unique 
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Table 3. NRE43-d8 motifs detected in nitrogen-related promoters from gene expression Cluster 1 



Maize annotation 


Description 




NRE43-d8 
motif 




Downstream nucleotides 


Percentage 
identity 


GRMZM2G010280 


High affinity nitrate transporter 
tnrtz. \} 


-319 


TGATCCTT 


-313 


GGCTGATCCC ACGGGATGAG GCCAAGCCCA 


93.24% 


-53 


CGACCCTT 


-47 


CATGTCCATG ACACGCCAGA GCTCAATCTT 


100.00% 


GRMZM2G453320 


Nitrate and chloride transporter 


-424 


CGACCTTT 


-418 


ATGATTTTGG GTCTTCTTTT TGAAAACGAA 


96.65% 


GRMZM2G080045 


Ammonium transporter 


-292 


CGACCCTT 


-286 


TTGCCGTGCG CTGCTCGAGT CTGCCTAACC 


100.00% 


-223 


CGTGCCTT 


-217 


CCGTTTTAGG TTTGATTCGT CGACTTGAAT 


90.03% 


-95 


CGACCATG 


-89 


TTCTTGTCAT CTCTTCAGAA CAGTCTGAAG 


91.12% 


GRMZM2G082343 


Symbiotic ammonium transporter 


-493 


CCACCGTT 


-487 


AGATCGGTCA CCAGGTCATA GTCCACCATG 


91.12% 


-303 


CGATCGTT 


-297 


CGTCGGGCGA TTGTTTATCC CCGGACTAAA 


95.61% 


unlVI/JVIzuUzo/ jo 


Ammonium transporter (0samt1.3) 


-208 


CGACACTT 


-202 


TATTGTAATT TTGGACTAGT CTCTCTTTTT 


94.29% 


-146 


CGATCTTT 


-140 


CTACAGTGCA AGATAATAAT GGAGTATCTC 


95.61% 


GRMZM2G175140 


Ammonium transporter (Osamtl.1) 


-169 


CGGCCCTG 


-163 


GGATCCTGGC CACCGTGGGT GGGCAGATTC 


90.67% 


-111 


CGATCCGT 


-105 


TGTTTGTTTT GCCGAATCAA AACTGCAATT 


93.24% 


GRMZM2G335218 


Ammonium transporter 
(Osamt2.2;Osamt3.1;Osamt5.1) 


-160 


CGATCCTT 


-154 


CTCTTCTCTC CTAGAGCCAC TCACCGGCGC 


98.95% 




Maize transporter consensus 




CGACCCTT 




Dicot- nitrite reductase consensus 
(Konishi and Yanagisawa, 2010) 




tGACcCTT 



probes analyzed in this study was 46 (out of 65) but we have 
presented data for all probes analyzed. 
Other study limitations 

On a cautionary note, the qRT-PCR data presented here 
(Fig. 2, Fig. S2) suggests that nitrogen-related gene expression 
in maize vegetative tissues may be more complex or nuanced 
than the 4 clusters identified by microarray analysis. For exam- 
ple, whereas Nar2.2 and a high affinity nitrate transporter 
(GRMZM2G455124) were grouped together into Cluster 4, the 
qRT-PCR expression data suggests overlapping but distinct pat- 
terns of expression: Nar2.2 peaked only in reproductive stage 
leaves (Fig. 2G) while the transporter peaked in both reproduc- 
tive stage leaves and roots (Fig. 2H). Another significant limi- 
tation of our study was that it excluded a significant number 
of nitrogen-related genes (see Methods for explanation) includ- 
ing genes encoding: nitrite reductase (which is induced in both 
maize leaves and roots); 44 mitochondrial glutamate dehydroge- 
nase (GDH); and 4 of the 6 known maize glutamine synthe- 
tases (GS) (GS2; GS1: Glnl-1; Glnl-2, Glnl-3). Therefore, 
the expression clusters described here should not be viewed as 
comprehensive. Finally, as the tissue sampling did not focus on 
senescing leaves, it likely caused probes related to scavenging and 
export of nitrogen to growing tissues (e.g., GS1), 45 to be under- 
represented in this study. 

Physiological relevance of each vegetative expression cluster 
In the context of the above study limitations, an intriguing 
observation from this study is that there appears to be a higher- 
order physiological pattern to the expression clusters (Fig. 1). 
Specifically, each of the 4 vegetative expression clusters appears 
to be selectively over-represented for genes responsible for spe- 
cific steps of the nitrogen uptake and assimilation pathway: 



Cluster 1: In this cluster, which showed peak expression in 
juvenile stage roots and leaves, 7/8 putatively unique probes 
corresponded to nitrate/nitrate-chloride transporters (4/8) and 
ammonium transporters (3/8), rather than to nitrogen assimi- 
lation genes (Fig. 1; Table S3). qRT-PCR data confirmed that 
expression of both a low affinity transporter (Nrtl.l) and a high 
affinity transporter (Nrt2.1) peaked in juvenile roots (Fig. 2). 
It was logical that both nitrate and ammonium transporters 
showed high expression, as the plants were fertilized with both 
sources of inorganic nitrogen throughout their life cycle. Given 
that we identified a distinct set of over-represented motifs in the 
promoters corresponding to Cluster 1 genes (Tables 2 and 3; 
see below), together these results suggest that maize may have 
a transcriptional module devoted to coordinating early stage, 
whole plant, inorganic nitrogen uptake and transport, possibly 
including transport of nitrogen into the vacuole (e.g., nitrate- 
chloride transporter), and loading/unloading of nitrate and 
ammonium to/from xylem. 40 The existence of an early root/ 
shoot inorganic nitrogen expression module is consistent with a 
recent study" which demonstrated that there is a peak in nitrate 
uptake and transporter expression in young maize roots when 
normalized to root dry weight, which results in an early stage 
peak in leaf nitrate concentration. Relative to root growth, the 
early stage peak was shown to correspond to a period of rapid 
shoot growth and hence increased nitrogen demand. 11 

Cluster 2: This cluster was leaf selective, with peak expres- 
sion generally at earlier developmental stages (Fig. 1). A closer 
examination showed that 9/11 putatively unique probes in 
Cluster 2 corresponded to proteins with nitrate transport activ- 
ity (5/11, including 1 unique nitrate-chloride transporter) as 
well as nitrate reductase activity (4/11) (Fig. 1; Table S3). qPCR 
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data confirmed that nitrate reductase gene Nr2 was leaf selective, 
with peak expression in VI stage leaves (Fig. 2C) . One interpre- 
tation of this data is that maize has an early stage leaf-selective 
module that is minimally responsible for unloading nitrate from 
xylem into leaf cell cytoplasm, followed by storage of nitrate in 
leaf vacuoles or conversion of nitrate to nitrite in the cytoplasm. 4 
Our data would suggest that these two steps are transcription- 
ally coordinated, perhaps facilitated by the motifs that were 
over-represented in the promoters associated with this module 
(Table 2). Maize leaf NR activity is associated primarily with 
mesophyll cells, with dramatic upregulation during greening. 46 

The expression of additional probes in Cluster 2, correspond- 
ing to 1 nitrite transporter and 1 ferredoxin-dependent GS (note: 
nitrite reductase was excluded from this study), suggests that 
other components of the nitrate assimilation pathway may also 
be part of this regulatory module. Given the existence of Cluster 
4, which was selective for leaf reproductive stages, an attrac- 
tive hypothesis is that Cluster 2 is coordinated with Cluster 1 
to satisfy the high shoot/root ratio during early development. 11 
As noted above, Cluster 1 was over-represented for probes asso- 
ciated with early stage root and leaf nitrate (and ammonium) 
uptake and transport, but not assimilation (Fig. 1), and hence 
the 2 clusters may be complementary. 

None of the Cluster 2 probes was associated with ammonium 
transport, even though 9 putatively unique probes correspond- 
ing to ammonium transporters were included in the study, and 
the plants were fed ammonium (Fig. 1). This result is consistent 
with previous observations demonstrating that de novo ammo- 
nium assimilation occurs primarily in the root. 40,47,48 

A final interesting observation for Cluster 2 is that the single 
ferredoxin-GOGAT affiliated probe in this study was expressed 
in this cluster (Fig. 1). This result is consistent with ferredoxin- 
GOGAT activity being dominant in leaf chloroplasts compared 
with NADH-GOGAT, which is expressed primarily in non- 
photosynthetic tissues. 17 ' 49 

Cluster 3: Cluster 3 was the root selective cluster that was 
expressed both early and late in development (Fig. 1). Cluster 3 
included 12 putatively unique probes, of which 3 probes corre- 
sponded to proteins with nitrate transporter activity and 1 with 
nitrate reductase (NR) activity. qRT-PCR results confirmed the 
root-selective expression of the high affinity co-transporter gene 
Nar2.1 in Cluster 3 (Fig. 2E; Table S5). In maize, though a con- 
siderable fraction of grain nitrogen originates from scavenged 
nitrogen from senescing leaves, 30—70% comes from post- 
silking uptake of soil nitrogen, especially nitrate. 45,50 In fact, a 
recent study by Garnett et al. 50 showed that there is a reproduc- 
tive (anthesis) phase peak in root nitrate uptake associated with 
increased expression of nitrate transporters. Enzymatic nitrate 
reductase activity has been measured in maize roots, with peak 
levels occurring in root tips. 51 

The remaining 8 probes in Cluster 3 corresponded to genes 
encoding the ammonia assimilation pathway (which was miss- 
ing in root Cluster 1), including: 1) 3 putatively unique probes 
corresponding to proteins with ammonium transport activity 
(these might enable soil uptake of ammonium, root vacuolar 
storage of ammonium, 52 or cytoplasmic to plastid transport to 



facilitate the first step of ammonium assimilation); 2) 1 or 2 
probe (s) corresponding to the next step in ammonium assimi- 
lation, specifically incorporation of ammonium with glutamate 
to form glutamine, catalyzed by glutamine synthetase (GS); 
3) 3 putatively unique probes corresponding to glutamate syn- 
thase (NADH-GOGAT); and 4) 1 amino acid transporter, 
potentially to transport glutamate, glutamine, or other amino 
acids. There were no probes expressed in Cluster 3 that corre- 
sponded to proteins with nitrite transport activity, while probes 
for nitrite reductase were not included in the study as already 
noted. 

Given that our study probes did not represent the complete 
set of nitrogen-related genes in the maize genome, one interpre- 
tation of Cluster 3 is that maize has a transcriptionally coor- 
dinated set of genes that acts in early-to-late stage maize roots 
and which encode the complete nitrate and ammonia uptake 
and assimilation pathway. As already noted, the existence of an 
ammonia assimilation pathway in roots is consistent with previ- 
ous studies demonstrating that the majority of de novo ammo- 
nium assimilation occurs in roots. 48,51,53 With respect to nitrate, 
assimilation of nitrate is reported to be more prevalent in leaves 
than in roots. 48 However, whereas feeding studies using maize 
roots showed a faster rate of amino acid accumulation when 
the nitrogen source was 15 NH 4 + rather than 15 N0 3 ", the differ- 
ence was thought to be caused by a lag period of induction of 
nitrate uptake rather than a slower steady-state rate of nitrate 
assimilation. 52 

As noted above, ferredoxin-GOGAT is reported to be pri- 
marily expressed in leaf chloroplasts, while NADH-GOGAT 
is expressed in non-photosynthetic tissues including root 
plastids. 17,40 Consistent with these observations, we observed 
the expression of 3 putatively unique probes corresponding 
to NADH-GOGAT genes in Cluster 3, but none for ferre- 
doxin-GOGAT (Fig. 1). qRT-PCR results confirmed that an 
NADH-GOGAT gene belonging to Cluster 3 was expressed in 
a root-selective manner (Fig. 2F; Table S5). 

Cluster 4: This late developmental stage, leaf selective clus- 
ter, which also included probes that were expressed in late stage 
roots (Fig. 1), was biased for probes corresponding to the nitrate 
transport and assimilation pathway (12/15 putatively unique 
probes). Specifically, Cluster 4 included: 1) 7 putatively unique 
probes corresponding to nitrate transporter activity, perhaps 
enabling nitrate unloading from xylem to leaf cytoplasm, or 
transport of nitrate into vacuoles for storage; 2) 3 putatively 
unique probes corresponding to nitrate reductase; and 3) 2 
probes corresponding to nitrite transporters. qRT-PCR con- 
firmed that a high affinity co-transporter (Nar2.2) was upregu- 
lated in reproductive stage leaves (Fig. 2; Table S5). Combined 
with the observation that Cluster 4 was over-represented for 
specific candidate promoter motifs (Table 2), these data suggest 
that maize has a transcriptional module that acts to facilitate 
reproductive stage nitrate transport and assimilation in leaves, 
perhaps to supply grain demand for nitrogen. 45,48,50 

With the exception of one NADH-dependent GS probe, 
it is curious that probes corresponding to downstream nitrate 
assimilation genes were not apparently expressed in Cluster 4; 



e26056-10 



Plant Signaling & Behavior 



Volume 8 Issue 10 



possible reasons include their exclusion from the experiment (see 
Methods) and/or that the downstream genes belong to a differ- 
ent regulatory module. In support of the latter hypothesis, nitrate 
reductase and nitrite reductase activities have been shown to be 
restricted to mesophyll cells in maize leaves, whereas glutamine 
synthetase (GS) and glutamate synthase (GOGAT) activities 
occur in both mesophyll and bundle sheath cells. 46 Furthermore, 
whereas the activities of nitrate reductase and nitrite reductase 
are upregulated during greening, greening causes minor effects 
on GS and GOGAT activities. 46 One interpretation of these 
results, when combined, is that the early (cytoplasmic) vs. late 
(chloroplast) steps in leaf nitrate assimilation belong to separate 
transcriptional modules, a hypothesis that will require further 
investigation. 

As already noted in the Study Limitations section above, a 
limit to our interpretation of Cluster 4 is that we under-sampled 
senescing leaves. In senescing leaves, ammonium is produced by 
protein degradation, and is subsequently assimilated into the 
amide group of glutamine by GS, for subsequent transport to 
growing tissues including grain. 54 Two ammonium transporters 
were expressed in Cluster 4, which could potentially be involved 
in the export of ammonium from the vacuole into the cytosol 
following proteolysis, as a substrate for GS1. 4 Another possible 
study limitation is that we sampled the second leaf above the 
ear to collect reproductive stage leaf data (Table SI). During 
maize grain filling, cytosolic GS1 isoforms have been shown 
to be expressed in the leaf below the ear, which is dedicated 
to ear feeding. 55 The single NADH-dependent glutamate syn- 
thase probe (ZmN-65; Table 1) observed in Cluster 4 leaves 
was highly expressed during the reproductive stages. The cor- 
responding gene could play an important role in nitrogen remo- 
bilization from leaf to grain. 

Promoter motifs underlying stage-selective gene expression 
A prediction of this study is that each vegetative nitrogen 
gene expression cluster is associated with distinct motifs that are 
over-represented in the associated set of promoters (Table 2). Of 
particular interest was the 8 bp motif (consensus CGACCCTT) 
that was over-represented in the promoters of Cluster 1 genes and 
which matched a sequence within the 43 bp nitrogen response 
element (NRE) (Table 3); the NRE was previously identi- 
fied in the promoter of the Arabidopsis nitrite reductase (NiR) 
gene, and shown to be both necessary and sufficient for nitrate- 
induced gene expression. 15 Within the 43 bp NRE, an alignment 
of monocot-dicot NiR promoters suggested that an 8 bp distal 
element [consensus (c/t)GaCcCTT] and a 5 bp proximal ele- 
ment [AAG(a/g)], separated by a 10 bp spacer, were somewhat 
conserved in the -100 bp to -240 bp promoter region, 20 with 
the distal 8 bp element being dominant. 19 Our predicted 8 bp 
Cluster 1 promoter motif matched the dominant distal element 
of the NRE and was similarly located in the -200 to -250 bp 
regions of Cluster 1 promoters (Tables 2 and 3). However, we 
failed to find the 5 bp proximal element in our study, as the 
Promzea software only searches for motifs > 6 bp. Recently, a 
yeast 1 hybrid screen using the 43 bp Arabidopsis NRE as bait 
identified the NIN-LIKE PROTEIN 6 (NLP6) transcription 
factor (and perhaps related NLP proteins) as being responsible 



for binding and nitrate-inducible activation of NRE in vitro and 
in planta. 21 Linker-scanner mutagenesis showed that the 8 bp 
distal element was especially responsible for nitrate activation by 
NLP6, in Arabidopsis protoplasts. 21 NRE-like sequences includ- 
ing the 8 bp element were also previously identified in the 3' 
region of the Arabidopsis nitrate reductase NIA1 gene 23 where 
it was shown to help mediate nitrate activation of this gene by 
NLP6. 21 In Arabidopsis, the NRE was only reported in nitrite 
reductase (NIR1) and nitrate reductase (NIA1) genes, 21,22 but 
in our study it was found to be over-represented in the promot- 
ers of genes encoding nitrate transporters (NRT1, NRT2.1) and 
ammonium transporters (Table 3), a novel observation that sug- 
gests that NLP6-like proteins may activate these promoters in 
maize. Here we refer to this conserved distal 8 bp motif as the 
NRE43-d8 motif. As NRE43-d8 was only over-represented in 
Cluster 1 promoters, one possibility is that it is an important 
cw-acting element of nitrogen-related genes during early-stage 
vegetative development in maize. 

Prior to the discovery of the NRE 43 bp motif noted above, 15 
an even more distal motif, also called the NRE [core A(c/g) 
TCA], was shown to be conserved in the promoters of nitrate 
reductase (NR) and NiR genes in both dicots and monocots 
including maize, but not tested for sufficiency of nitrate induc- 
tion. 56 In maize, several of the distal NREs [A(g/c)TCA] were 
located in the -500 to -1000 bp promoter region, 56 more distal 
than the region analyzed in this study, perhaps explaining why 
we did not recover this motif in our search. 

Our analysis did however reveal another putative nitrogen- 
related motif of interest, the G-box 10 element, which was 
over-represented in Cluster 1 gene promoters (Table 2). The 
core of the G-box 10 element is identical to the E-box motif 
(CANNTG) which is necessary for nitrate induction of the 
NIA1 (nitrate reductase) promoter in Arabidopsis. 24 The E-box 
motif and other non-NRE motifs that were identified in this 
study (Table 2) may be the yet to be identified binding sites for 
NLP7, a regulatory factor recently shown to bind 851 genes in 
Arabidopsis within 10 min of nitrate addition including genes 
encoding nitrate reductase (NIA1), GS2, and transporters 
NRTl.l, NRT2.1, and NAR2. 22 Interestingly, nitrate triggers 
rapid nuclear retention of NLP7 leading to chromatin binding. 22 
NLP6 and NLP7 are related, but they have been reported to 
have different chromatin binding sites. 22 

Significance for crop improvement 

It is estimated that only 30-80% of nitrogen fertilizer is 
taken up by maize roots, with the remainder leached or vola- 
tilized. 3,57 " 5 ' However, improving the nitrogen use efficiency of 
maize will require a better understanding of nitrogen uptake 
and assimilation in different tissues and stages of development. 
For example, in newer maize hybrids, an increase in nitrogen 
uptake efficiency has been attributed specifically to increased 
post-flowering nitrogen uptake. 60 Unfortunately, the underly- 
ing transcriptional regulation of nitrogen-related genes dur- 
ing maize development has been poorly understood, with the 
first comprehensive study of nitrate transporters only recently 
reported. 11 Here we identified distinct developmental modules 
of nitrogen-related gene expression in maize. We have discussed 
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the possibility that each module may be fine-tuned for the dif- 
ferent nitrogen uptake and assimilation needs of each tissue and 
stage of development. Clustering of nitrogen-related genes by 
their developmental tissue and stage of expression appears to 
have assisted in the identification of putative cw-acting promoter 
motifs; this strategy might represent a path forward in other 
plants. These advances may assist in targeting long-term efforts 
for improving nitrogen use efficiency in maize. 

Materials and Methods 

Plant growth and tissue harvest 

Syngenta inbred SRG200 seeds were grown in a greenhouse 
during the summer of 2007 at the University of Guelph, using 
the following conditions: 16 h light (-600 ixmol m" 2 s" 1 ) at 28°C, 
8 h dark at 23°C, and 50% relative humidity. Plants were grown 
semi-hydroponically in pots containing Turface® clay, and irri- 
gated with a nutrient solution containing: 0.4 g/L 28-14-14 fertil- 
izer [28% total N (1.6% nitrate, 0.4% ammonium, 26% urea), 
14% P 2 0 5 , 14% K 2 0], 0.4 g/L 15-15-30 fertilizer [15% total N 
(nitrate 8.8%, ammonium 2.95%, urea 3.25%), 15% P 2 0 5 , 30% 
K 2 0], 0.2 g/L NH 4 N0 3 , 0.4 g/L of MgS0 4 «7H 2 0, and 0.03 g/L 
of micronutrient mix (S, Co, Cu, Fe, Mn, Mo and Zn). Three bio- 
logical replicates per tissue/stage were harvested, always at -11 AM. 

Microarray analysis 

Details of the microarray and hybridization procedure have 
been described previously. 27 Briefly, RNA was isolated from 
50 tissues/stages (Table SI) and hybridized onto maize B73 
Affymetrix 82K Unigene arrays. The Syngenta GeneChip micro- 
array was custom-designed in 2005. The starting array consisted 
of approximately 87 000 probe sets representing 82 000 Zea mays 
unigenes (as defined by Pubmed) or expressed sequence tag (EST) 
clusters, which were trimmed for further analysis (Table S2), 
resulting in only a subset of nitrogen-related genes to be reported 
in this study (see details below). Each probe set consisted of an 
average of 16 perfect-match probes, of 25 mer oligonucleotides. 
For each tissue-stage, 3 hybridizations were conducted using 
RNA isolated from 3 biological replicates representing 3 pools of 
plants. Total RNA was isolated and used for cDNA, cRNA syn- 
thesis, and labeling using standard Affymetrix protocols. Labeled 
cRNAs were then fragmented and applied to the maize custom 
GeneChip microarray for molecular hybridization. 61 The array 
images with hybridization signals were acquired and quantified 
by GeneChip Operation System (GCOS) software (Affymetrix). 
The quality of the hybridizations was assayed using Expressionist 
(GeneData). Experiments were repeated for any arrays that ini- 
tially failed to pass the quality assays. Array expression was nor- 
malized using the RMA method 42 from Bioconductor. 62 The 
significance of differential expression for the 82 K probes was 
measured by ANOVA using R. The P values of relative microar- 
ray expression were adjusted using the Bonferroni method. 63 
Annotation and clustering of nitrogen-related genes 
Clustering was conducted using K-means clustering 64 of the 
relative expression data after RMA normalization. The number K 



of clusters was estimated by testing the uniformity of the cluster 
for K from 2 to 50 (for comparison of 50 tissues) or 17 (for com- 
parison of leaf and root stages only). Array probes corresponding 
to nitrogen-related genes were retrieved using 3 methods. First, 
as the probes were designed from the maize Unigene set, the cor- 
responding original GenBank sequences were used to retrieve 
matches, using nucleotide BLAST searches against the B73 maize 
genome (MaizeSequence.org, release 4a. 53). Probe sets with no 
expression (relative expression < 100) in any of the 150 microar- 
ray experiments were removed from the annotation (26 989 probe 
sets). If 75% of the probes in the probe set (12/16) matched the 
same gene model, the probe set was identified as a match for that 
gene. If only 12%-< 75% of the probes in the probe set (1/16- 
11/16) matched the same gene model, the probes were considered 
as a partial match to that gene. A probe was required to have 85% 
sequence identity with the gene model to be considered as a valid 
match of that gene. A total of 33 664 probe sets that matched to 
unique gene models were mapped using these steps (Table S2). 
Exonerate alignment 65 was used to annotate an additional 9919 
probes. Nucleotide BLAST was not successful in identifying 
EST matches because of the biases created by gene model issues 
related to gene-calling software (GeneBuilder or FGENESH). 
The remaining 12089 probe sets on the array showed expres- 
sion, but did not map to the maize genome. After the elimination 
of the probe sets with either no expression, cross-hybridizing, or 
redundant probe sets, there were 22787 high-quality annotated 
probes. Subsequently, the array probe sequences were re-screened 
for matches with EST sequences from NCBI using BLAST. 
Finally, nitrogen uptake and assimilation keywords were used to 
search gene annotations and protein domains in the B73 maize 
genome from MaizeSequence.org, and the search results were 
each matched to a GenBank protein with the greatest homology 
with the microarray probe set. After this pipeline, probes corre- 
sponding to nitrogen-related genes such as nitrite reductase were 
found to be missing which is a limit of the current study. 

Quantitative real-time reverse transcription PCR 
(qRT-PCR) 

The genes and PCR primers used for qRT-PCR validation are 
listed (Table S4). An aliquot of the RNA used for microarray 
hybridization was quality control tested using the Agilent 2100 
Bioanalyzer. Total RNA was used to create cDNA. Subsequently, 
qRT-PCR was performed using Perfecta SYBR® Green FastMix 
ROX™ (Quanta BioSciences, Inc, Gaithersburg, MD USA) and 
a StepOnePlus™ Real-Time PCR System (Applied Biosystems- 
Life Technologies Corp, Carlsbad, CA USA). Amplification 
conditions were 95°C for 3 min, followed by 40 cycles of: dena- 
turation, 95°C for 15 s; annealing (55°C for NAR2.1 and 60°C 
for the other primers) for 30 s; extension at 72°C for 1 min. The 
results were normalized using the Livak method to report relative 
expression. 66 

Promoter motif discovery 

Microarray probes were matched to known genes from the 
filtered B73 maize gene set from MaizeSequence.org as previ- 
ously described. 28 Over-represented de novo motifs that were 
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present in the promoters (-500 to +1 region) of the genes 
within each expression cluster were then searched using a Perl- 
based motif discovery program developed by us for the maize 
genome, called Promzea. 28 Briefly, Promzea filters and com- 
bines results from 3 motif discovery tools (Weeder, 67 MEME, 68 
and BioProspector 69 ) and performs additional validation tests. 
Retrieved motifs are subsequently matched to previously 
defined motifs using Athamap 70 and PLACE 71 using STAMP 
software. 72 Promzea predicts motifs > 6 bp in length. 
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